Mon. Not. R. Astron. Soc. 000, [T]-?? (2009) Printed 5 August 2009 (MN style file v2.2) 



Circumstellar interaction in type Ibn supernovae and 
SN 2006jc 



O N.N. Chugai * 

Institute of Astronomy, RAS, Pyatnitskaya 48, 119017 Moscow, Russia 

^ : 

• Accepted 2009 Received 2009; in original form 2009 

in 



O 

u 

o 



> 
00 

in 
o 



00 

o 
o 



X 



ABSTRACT 

I analyse peculiar properties of light curve and continua of enigmatic Ibn super- 
novae, including SN 2006jc, and argue in favour of the early strong circumstellar inter- 
action. This interaction explains the high luminosity and fast flux rise of SN 1999cq, 
while the cool dense shell formed in shocked ejecta can explain the smooth early con- 
tinuum of SN 2000er and unusual blue continuum of SN 2006jc. The dust is shown 
to condense in the cool dense shell at about day 50. Monte Carlo modelling of the 
He I 7065 A line profile affected by the dust occultation supports a picture, in which 
the dust resides in the fragmented cool dense shell, whereas He I lines originate from 
circumstellar clouds shocked and fragmented in the forward shock wave. 
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• 1 INTRODUCTION 



Recently a new family of supernovae (SN), so called SN Ibn, 
came in the limelight (Matheson et al. 2000; Pastorello et 
al. 2008). These objects, showing strong narrow (FWHM~ 
2000 km s~^) He I emission lines and week or no Hq emis- 
sion, are suggested to be SN Ibc interacting with a helium- 
rich circumstellar matter (CSM) (Foley et al. 2007; Pas- 
torello et al. 2007; Pastorello et al. 2008). Four known mem- 
bers of this family are: SN 1999cq (Matheson et al. 2000), 
SN 2000er and SN 2002ao (Pastorello et al. 2008), and well 
studied SN 2006jc (Nakano et al. 2006; Foley at al. 2007; 
Pastorello et al. 2007; Di Carlo et al. 2008; Smith et al. 
2008; Anupama et al. 2008). In SN 2006jc the CS inter- 
action is manifested also by X-ray emission (Immler et al. 
2008). At the position of SN 2006jc an optical transient was 
detected in Oct. 2004 with the maximal absolute magnitude 
of Mr ~ -14 (Nakano et al. 2006; Foley et al. 2007; Pas- 
torello et al. 2007). This outburst is blaimed for the mass 
ejection by presupernova (Foley at al. 2007; Pastorello et al. 
2007). 

A striking feature of SN 2006jc is the early dust for- 
mation {t ~ 50 d) indicated by the infrared emission of 
hot (T ~ 1800 K) dust after about day 50 (Di Carlo et al. 
2008; Smith et al. 2008; Mattila et al. 2008) and blue shift 
of emission hues (Smith et al. 2008; Mattila et al. 2008). 
Although dust might form in the unshocked ejecta as shown 
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by Nozawa et al. (2008), the dust formation in the shocked 
cool dense shell (CDS) seems to be preferred (Mattila et 
al. 2008). In fact, dust condensation in the CDS has been 
invoked already for type IIL supernova 1998S (Pozzo et al. 
2004). 

A general wisdom is that the light of SN 2006jc is pow- 
ered by the radioactive decay of ^'^Ni similar to normal 
SNe Ibc (Tominaga et al. 2008). The X-ray emission is ex- 
plained by Tominaga et al. (2008) as a result of interaction of 
SN Ibc ejecta with smoothly distributed rarefied CSM with 
the total mass of ~ 0.003 Mq. In this model both reverse 
and forward shocks are adiabatic, i.e., the dust formation in 
the shocked gas is precluded. In scenario proposed by Mat- 
tila et al. (2008) SN 2006jc ejecta collide with a massive 
~ 0.5 Mq CS shell at the distance of 1.5 x lO^'^ cm. The 
massive CDS formed in the forward shock was presumably 
the site where the dust condensed after day 50 (Mattila et 
al. 2008). Unfortunately authors did not discuss a possible 
contribution of CS interaction in the bolometric luminosity. 

The present study is motivated by the unsettled issue of 
the dust-forming site and striking facts presently escaped the 
detailed analysis, viz., (a) luminous maximum and fast pre- 
maximum flux rise of SN 1999cq; (b) smooth early contin- 
uum of SN 2000er; and (c) odd blue continuum of SN 2006jc. 
The goal of the paper is to propose a model that could ac- 
count for the above properties and some other observational 
data, including the dust formation. The paper is organized 
as follows. Arguments in favour of a strong CS interaction in 
SNe Ibn at the early stage are considered in Section (2] The 
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interaction models with and without ^^Ni are then applied 
to account for the bolometric light curve, X-ray emission, 
and moderate expansion velocity of outer layers (Section |3]). 
The parameters of the CDS recovered in these models are 
used to study the dust formation in terms of the theory of 
homogeneous nucleation (Section The Monte-Carlo sim- 
ulations of dust occultation effects in He I lines are made to 
illustrate that the proposed general picture is sensible (Sec- 
tion The results and some worrisome issues are discussed 
in Section [6l 

Hereafter we adopt for SN 2006jc the explosion date 
2006 Sep. 21 or JD=2454000 (Mattila et al. 2008) and the 
distance of 26 Mpc according to redshift from the Lyon Ex- 
tragalactic Database and assuming Ho = 70 km s"^ Mpc"\ 
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Figure 1. Blue continuum of SN 2006jc on day 22. The model 
of line emission from cool dense shell (thick solid line) with the 
reddening A{V) = 0.15 is overplotted on the observed spectrum 
(Fooley et al. 2008). 



2 EARLY CS INTERACTION 

Pastorello et al. (2008) emphasise the remarkable homogene- 
ity in the spectral properties among SNe Ibn. This fact pro- 
vides us a confidence that the conclusions made from the 
analysis of particular objects can be applied to other mem- 
bers of the family. Here I consider observational data which 
suggest a dominant contribution of CS interaction in the 
early luminosity and spectra of SNe Ibn. 



2.1 Light curve and CSM 

Unusually high luminosity of SN 1999cq (type Ibn) at the 
light maximum with the absolute magnitude Me. ~ —19.9 
mag (Matheson et al. 2000; PastoreUo et al. 2008) and the 
fast pre-maximum flux rise by > 3 mag in four days (Math- 
eson et al. 2000) are unprecedented for SNe Ibc including 
luminous SN 1998bw. Both high luminosity and fast bright- 
ening suggest that early SN 1999cq was powered by the 
ejecta interaction with a dense CSM, a scenario similar to 
that proposed for the early light curve of type IIL SN 1998S 
(Chugai 2001). SN 2006jc was detected on day 17 after the 
adopted explosion date (Nakano et al. 2006; Pastorello et al. 
2008) and its absolute magnitude at the discovery was simi- 
lar to that of SN 1999cq at the same age. It seems plausible, 
therefore, that the prediscovery light curve of SN 2006jc was 
similar and the early luminosity of SN 2006jc was also af- 
fected by the CS interaction. The early interaction should 
lead to the deceleration of outer layers of ejecta which is 
supported by the fact that the early {t < 10 d) He I 5876 
A emission of SN 2000er shows the maximal expansion ve- 
locity of only ~ 9000 km s"^ (Pastorello et al. 2008), rather 
moderate value for early SN Ibc. 

The total energy radiated by SN 1999cq during the first 
10 days after the light maximum is about Er ^ 2 x lO'*^ erg 
according to the available light curve (Matheson et al. 2000; 
Pastorello et al. 2008) . To produce that energy a shock wave 
with velocity of 10* km s~^ had to interact with the CS mass 
of about Ales ^ 0.02 Mq assuming 100% radiation efficiency. 
For that CSM mass enclosed within r ^ vt ^ 10^^ cm the 
density parameter for the steady outflow turns out to be 
w — Mcs/r ~ 4 X 10^'' g cm~^. For the density distribution 
in outer SN ejecta p oc v~'^ the shocked ejecta should be by 
2.5 times more massive (Chevalier 1982b), so the expected 
CDS mass on day 10 could attain ^ 0.05 Mq. 



2.2 Early smooth continuum of SN 2000er 

The earliest spectrum of SN 2000er shows a smooth contin- 
uum (Pastorello et al. 2008), very unlike bumpy continua 
of any early SN Ibc. I suggest that the smooth continuum 
forms in a thin CDS likewise in early SN 1998S (Chugai 
2001). The smooth continuum of SN 2000er grows bumpy 
on the time scale of flve days after the first spectrum (Pas- 
torello et al. 2008), which means that the continuum forms 
in a partially transparent regime, i.e., the optical depth is 
r < 1. For the luminosity of ~ 10"*^ erg s~^ on day 10 (Pa- 
storello et al. 2008) and photospheric radius r ~ ~ 10^^ 
cm the effective temperature of the continuum turns out to 
be ~ 10'' K. 

A question arises, whether a helium-rich CDS with a 
mass of ~ 0.05 Mq is able to maintain a smooth contin- 
uum at the age of ~ 10 days? To illustrate the possibility 
I consider the radiation transfer in the isothermal helium 
slab with Thomson scattering and absorption taken into ac- 
count. To facilitate estimates I adopt the LTE approxima- 
tion. With the boundary condition _F(0) = c[/(0)/V3 (U is 
the radiation density) is then 

47r 7^(1 - E) 

F.-— B.(r)^^^^^^=^^, (1) 

where E = exp(— r-^Se), t is the extinction optical depth 
and e is the thermalization parameter, i.e., absorption to 
extinction ratio; both r and e are functions of frequency. 
This expression permits us to calculate luminosity of the 
CDS for a certain gas temperature using Kramers-Unzold 
approximation for the absorption coefficient of the helium- 
rich matter. 

The density in the CDS is found from the pressure equi- 
librium adopting the wind density parameter w = Ax 10^^ g 
cm~^, shell velocity of lO** km s^^, and the radius of ~ 10^^ 
cm. To produce the continuum luminosity of ~ 10'''^ erg 
s~^ the equation ^ requires the gas of ~ 12500 K for 
the CDS mass of 0.05 Mq. In the range A > 4000 A the 
model spectrum is close to black-body with the tempera- 
ture T ~ 10200 K. Although the extinction optical depth of 
the shell at 6000 A is only ~ 0.4 the equilibration is rather 
efficient because thermalization parameter is large e > 0.92 
in the range A > 3000 A. We thus conclude that the smooth 
continuum may arise in early SN Ibn. The epoch of the early 



Circumstellar interaction in type Ibn supernovae and SN 2006jc 3 



smooth continuum of SN 2006jc preceeded the discovery and 
therefore was missing. 

2.3 Blue bumpy continuum at t > 20 d 

The smooth continuum of SN 2000er grew bumpy on a time 
scale of about one week (Pastorello et al. 2008), although 
the short period of SN 2000er observations does not per- 
mit one to follow the early evolution from smooth till fully 
developed bumpy continuum we see in SN 2006jc on day 
22. No doubt, the bumpy continuum is composed by nu- 
merous lines, presumably of Fell as suggested by Foley et 
al. (2007). However, unlike any normal SN Ibc, including 
SN 1998bw, the spectum of SN 2006jc is unusually blue: the 
colour B — R ^ —0.45 (Foley et al. 2007) compared with 
0.48 of SN 1998bw at the light maximum (Galama et al. 
1998). The spectrum also does not show broad absorption 
lines characteristic of SN Ibc. This indicates that the config- 
uration of "photospere" surrounded by an extended free ex- 
pansion atmosphere is irrelevant in this case. Following the 
picture invoked for the early smooth continuum I assume 
that the bumpy blue continuum of SN 2006jc arises from 
the same CDS but this time with the dominant emission in 
metal lines. This picture is similar to SN 2002ic (Chugai et 
al. 2004a), so to demonstrate the possibility I use the same 
model and the same line list of ~ 18000 lines of iron-peak 
elements in the range of 3500-10000 A. 

The calculated continuum and observed SN 2006jc spec- 
trum on day 22 (Foley at al. 2007) are shown in Fig. [1] 
The model suggests the emission of randomly oriented plane 
sheets distributed in the shell AR/R = 0.1. The total mass 
of cold fragments is 0.05 Mq and the solar mass fractions 
of iron peak elements (Ti, V, Cr, Fe, Ni) are assumed. The 
adopted expansion velocity is 9000 km s~^, ionization tem- 
perature 7000 K, excitation temeperature 8000 K, and the 
area ratio A = S/A-kR? — 20. Just to recall, the area ra- 
tio is the ratio of cumulative surface area of fragments to 
the surface area of the CDS; this parameter characterizes a 
mixing degree. The model spectrum is sensitive to the varia- 
tion of the excitation temeperature within 10% and in lesser 
extent to to the shell mass and thickness, which may vary 
within 30% without pronounced change of the spectrum. 
The adopted value of the shell thickness 0.1-R corresponds 
approximately to the width of mixing layer formed as a re- 
sult of a Ray leigh- Taylor instability of the thin shell (Cheva- 
lier & Blondin 1995). The model qualitatively reproduces the 
observed spectrum and bolometric luminosity ~ 3 x 10*^ erg 
s"^ (cf. Mattila et al. 2008). The high relative intensity of 
the blue continuum is the result of a larger number of lines 
per unit wavelength interval in the blue compared to the 
red part of the spectrum and a weak saturation of intensity 
in the blue for a moderate area ratio. Given satisfactory fit 
and sensible values of model parameters, the proposed pic- 
ture for the blue continuum provides us indirect support for 
the conjecture of the strong early CS interaction in SNe Ibn. 



2.4 CS interaction and emission lines 

Spectra of SN 2006jc (Foley at al. 2007; Pastorello et al. 
2007) suggest that the velocity of He I line-emitting gas lies 
in the range of 800 - 6000 km s"^ with the velocity of the 
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Figure 2. Emission lines in the red part of the spectrum of 
SN 2006jc on day 22. Thick line shows synthetic spectrum; thin 
lines is the observed spectrum (Fooley et al. 2008); dashed line 
is adopted continuum. Inset shows profiles used in the synthetic 
spectrum for Hel lines (dotted line) and metal lines (solid line) . 



bulk of emitting material in the range of 800 — 3000 km s~ . 
The minimal velocity of ~ 800 km s~^ could be associated 
with the cloud shocks speed and perhaps with the velocity 
of the undisturbed CSM. The latter is about 620 km s~^ as 
indicated by the absorption minimum of narrow lines of O I 
7773 A and Hel 6678 A (Anupamaet al. 2008). Interestingly, 
the expansion velocity of CSM of SN 2000er indicated by 
narrow Hel absorption components in early spectra is 800- 
900 km s~^ (Pastorello et al. 2008), close to this value. 

To account for the velocity spectrum of the Hel line- 
emitting gas I suggest that the bulk of the Hel emission 
originates from shocked CS clouds. Two-dimensional hydro- 
dynamic simulations of interstellar cloud engulfed by the 
blast wave (Klein et al. 1994) prompt us a picture in which 
the velocities of the line-emitting gas are produced by slow 
(because of density contrast) radiative shocks of CS clouds 
in combination with the fragmentation cascade of shocked 
clouds and acceleration of fragments. This scenario is similar 
to that proposed for the Ha emission in SN 2002ic (Chugai et 
al. 2004a) . The total optical luminosity of He I lines is about 
L ~ 5 X lO''" erg s'^ on day 22 (Smith et al. 2008). The cor- 
responding radiation efficiency of the forward shock kinetic 
luminosity is 77 = L/(0.5u;i'J). With the forward shock veloc- 
ity Vs — 9000 km s~^, wind density parameter w = 4,x 10^^ 
g cm~^ the required efficiency should be 77 ~ 0.003, a rea- 
sonable value. 

The other emission lines seem to originate from the 
same line-emitting gas. This is demonstarted by the syn- 
thetic spectrum of SN 2006jc on day 22 in which profiles and 
intensities were adjusted to fit the observed spectrum (Fig. 
[21 . Apart from Hel, OI, and Call lines the synthetic spec- 
trum includes Mgll lines (Pastorello et al. 2007; Anupama 
et al. 2008) and two C I lines 9088 A (^P°-^P) and 9406 A 
(^P°-^D), which are also pronounced in novae (Williams et 
al. 1991). There is a vague hint of a week N I 7452 A emis- 
sion with the intensity ratio NI 7452 A/O I 7773 A of ~ 0.1 
(Fig. [2| . The fact that profiles of all the lines are identical 
confirms that these emissions originate from the same line- 
emitting gas, which I associate with the shocked CS clouds. 
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Figure 3. Density distribution of the circumstellar gas for stud- 
ied models. Thin, intermediate and thick solid lines correspond 
to the model Al, A3, and A2 respectively, while the dotted line 
stands for the model B. 



3 INTERACTION MODEL 

The CS interaction and bolometric light curve are simulated 
using the model described earlier (Chugai 2001). To recap, 
the hydrodynamics of the interaction is treated in the thin 
shell approximation (Chevalier 1982a). The X-ray emissions 
of the reverse and forward shocks are calculated assuming a 
homogeneity of each post-shock layer and adopting the cool- 
ing function of Sutherland & Dopita (1993). The electron 
temperature in the postshock gas is calculated taking into 
account partial Coulomb equilibration. The absorbed and 
emergent X-ray luminosity is computed assuming the unab- 
sorbed emissivity spectrum to be e(-E) oc E~^'^ exp{—E/kT) 
and taking into acount the absorption by the cool gas (SN 
ejecta, CDS and unshocked CSM). The model despite many 
simplifying assumptions seems to be reasonable which is con- 
firmed by the modelling of SN 1999em X-ray flux (Chugai 
et al. 2007). The absorbed X-ray luminosity is reprocessed 
into the ultraviolet (UV), optical, and infrared photons. We 
adopt that the bulk of this radiation falls in the band in 
which the observed bolometric light curve is sampled. This 
assumption is sensible at the early epoch (t < 1 month), 
when the cool gas is opaque in UV region and rather crude 
at the late epoch (t ~ 1 yr), when significant fraction of the 
reprocessed radiation can escape in UV lines (Chevalier & 
Fransson 1994). The interaction luminosity combined with 
the diffusion luminosity of the ejecta produces the model 
bolometric luminosity. In turn, the diffusion luminosity is 
modelled using the Arnett (1980) approximation. 

A smooth CSM is assumed, neglecting dumpiness, 
which implies that the calculated luminosity and temper- 
ature of the X-ray emission from the forward shock perhaps 
are not reliable. I present two versions of the interaction 
model: with ^^Ni and without ^''Ni. The gamma-ray depo- 
sition is calculated assuming that ^^Ni is mixed in the inner 
90% of the ejecta mass. The ejecta density distribution is ap- 
proximated by the exponential law p oc exp{—v/vo), while 
the radial density distribution of the CSM is adjusted to de- 
scribe the bolometric and X-ray luminosities and the radii 
of the dusty shell associated with the CDS. The composi- 
tion of major elements is assumed to be He : C : O = 0.9 : 
0.06 : 0.04. This choice, called a "standard abundance", is 
similar to that of Mattila et al. (2008). The pre-SN radius 



-Ro = 10 Rq is assumed; the exact value is irrelevant, unless 
the radius is very large, _Ro 3> 10 Rq. 

The Table 1 shows model parameters: the explosion en- 
ergy, ejecta mass, ^''Ni mass, mass of the CS envelope, mass 
of the CDS formed in the reverse shock by day 50, and the 
wind density parameter w in the inner region r < 2 x 10^^ 
cm. Our models admit formation of the CDS with a mass of 
0.01 — 0.09 Mq in the forward shock during initial phase of 
5-15 days depending on the model. In a more realistic case 
of clumpy CSM the situation with the forward shock CDS 
would become more complicated. Radiative cloud shocks 
would result in the formation of cool CS gas in a forward 
shock. However, the expected kinematics and distribution 
of this cool gas would be essentially different compared to 
the model of a smooth wind. Models with ®®Ni differ by the 
ejecta parameters and the wind density. It should be stressed 
that ejecta mass and energy are poorly constrained by the 
light curve and spectrum. The high energy model A2 cor- 
responds to the preferred model of Tominaga et al. (2008), 
while the models Al and A3 have moderate explosion en- 
ergy. The CSM density distributions for all the models are 
presented in Fig. [3] 

All the models satisfactorily reproduce the observed 
bolometric light curve (Fig. 3^) including the maximal lu- 
minosity of SN 1999cq estimated from Mr magnitude (Pas- 
torello et al. 2008) . Note, we do not plot the low-mass model 
A3 since its behavior repeats that of the model Al in all the 
panels of Fig. ID Models with ^''Ni (Al and A3) reproduce 
the X-ray flux within la errors (Fig. HJj), while the model 
A2 with the high energy lies beyond la range. The model 
B predicts too high X-ray luminosity that exceeds by a fac- 
tor of ~ 5 X 10^ the observed X-ray luminosity around day 
50. This is a serious problem for the model without ^*Ni. 
Yet one cannot rule out that the dumpiness of CSM, ig- 
nored here, might modify the interaction zone in such a way 
that X-rays were efficiently absorbed. The observed X-ray 
fiux seem to show a drop by a factor 2.5 between the maxi- 
mum on day ~ 110 and the epoch of ~ 135 d (Immler et al. 
2008). This drop cannot be reproduced in any of our interac- 
tion model, although models Al and A3 are consistent with 
observations within la errors. If the drop is real then shocks 
in CS clouds engulfed by the forward shock layer contribute 
significantly to the X-ray emission. In this case the drop of 
X-ray flux could be related with the disappearence of CS 
clouds beyond the radius r ~ 10^^ cm. 

A massive (0.2 - 0.3 Mq) CDS forms by day 50 in aU 
the models (Table 1); the expansion velocity at the age of 
50-100 days is ~ 10^ km s"'^ in models Al, A3 and B and 
~ 26000 km s"^ in the high-energy model A2 (Fig. |4};). Ac- 
cordingly, for models Al, A3 and B the CDS radii are close 
to photometric radii of the dusty shell (Fig. UJl) reported 
by Mattila et al. (2008), while the high-energy model A2 
has essentially larger CDS radius. Interestingly, for the con- 
sidered models the CDS mass at t = 5 d falls in the range 
0.05 — 0.06 Mq, the value used in Section 2 to analyse the 
smooth continuum of SN 2000er. 

1 explored also the model in which the CSM with the 
mass of 0.1 — 0.5 Mq has a structure of a dense thin shell 
with the radius of 4 x lO'^^ cm. This choice corresponds to a 
scenario in which the CS shell was violently ejected with the 
velocity of 600 km s~^ during the outburst two years prior 
to the SN explosion. In this case the bolometric light curve 
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Table 1. Parameters of models of SN ejecta and CSM 



Model 
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M 






Mcds 


w 




10"'i erg 






Mq 




lO^S g/cm 


Al 


1.5 
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0.12 
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0.21 
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10 
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Figure 4. Left panels: Bolomctric light curves (a) and X-ray lu- 
minosity ib) escaping after absorption in unshocked ejecta, CDS 
and CSM for the models Al (thin solid line), A2 (thick solid line) 
and B (dotted line). Crosses show the observed bolometric light 
curve according to Pastorello et al. (2008); circle represents the lu- 
minosity of SN 1999cq at the light maximum; squares show X-ray 
luminosity in 0.2-10 keV band (Immler et al. 2008). Right panels: 
Velocity (c) and radius (d) of the cool dense shell. Triangles show 
radii of the dusty shell reported by Mattila et al. (2008). 



is notably affected by the CS interaction and the model light 
curve with or without ^^Ni is inconsistent with observations. 

To summarize, the bolometric light curve of SN 2006jc 
can be explained by models with and without ^®Ni; both 
cases suggest the dense wind at the distance r < 2 x 10^^ cm. 
In all the models massive CDS forms which might become 
the dust formation site. 



4 DUST FORMATION 

For the dust to condense the gas and dust temperatures 
should become lower than the condensation temperature 
(Tc). Unfortunately, the CDS gas temeperature (Tg) cannot 
be reliably determined from the radiative energy balance be- 
cause of complexities of molecular compounds and processes 
involved. We therefore adopt a sensible assumption that Tg 
is equal to the radiation temperature Tr — (F/ac)^^'^, where 
F is the radiation flux at the CDS radius, a is the radi- 
ation constant, c is the speed of light. This approximation 
probably is not perfect because the gas irradiated by diluted 
hot radiation is usually warmer than at the radiation tem- 
perature. To take into account a possible deviation from the 
radiation temeperature I also consider the case, in which the 
gas temperature is equal to the dust temperature computed 
in a formal way; this approximation presumably crudely re- 
flects the physics of heating and cooling of molecular gas. 
The dust temperature Td is calculated from the balance be- 



tween the absorption of SN radiation and the dust emission. 
We approximate the SN radiation by the diluted black body 
spectrum with the temperature of 10* K; dust grains are as- 
sumed to be graphite spheres with the radius of 0.03 fim and 
the absorption efficiency is set according to Draine & Lee 
(1984). Calculations show that on day 50 the dust tempera- 
ture turns out ~ 10% larger than the radiation temperature. 

Civen Tg < Tc, the dust, nevertheless, forms only, if 
a characteristic time of the dust formation is essentially 
smaller than the expansion time. The time scale of the dust 
formation is determined by the nucleation rate and growth 
of nuclei due to the gas accretion onto the grains. The rate of 
the dust growth can be easily estimated using results of the 
interaction modelling. At the age t — 50 d the carbon num- 
ber density in the CDS for the explored models with stan- 
dard abundance is n ~ (2 — 10) x 10^^ cm~^ assuming that 
all the oxygen is bound in CO. Adopting Tg — 1800 K and 
sticking parameter of 0.5 we find that it takes ~ (1—3) x 10^ s 
for the dust grains to grow a ~ 10~^ cm. This shows that 
the duration of the dust condensation probably is not con- 
strained by the dust growth. 

In line with the theory of homogeneous nucleation (cf. 
Drain & Salpeter 1977; Hasegawa & Kozasa 1988), the nu- 
cleation time scale can be expressed via the number density 



where n„ 



(2) 



nn/(47raV3) and ^ ~ 2 x 10" 



volume per carbon atom in a condensed phase. The nucle- 
ation rate (Draine & Salpeter 1977; Hasegawa & Kozasa 
1988) is defined by the accreation rate onto clusters close to 
equilibrium point 



J: 



a^:'Kr\uZr? exp(- 



167va /3g kT) , 



(3) 



where a is the sticking parameter (we adopt a = 0.5), Zel- 
dovich factor Z takes into account non- equilibrium distri- 
bution of nuclei sizes, u is the average thermal velocity of 
carbon atoms, Vc = —2a/g is the radius of a cluster of criti- 
cal size, a — 1400 erg cm~^ is the energy of graphite surface 
tension (Tabak et al. 1975), g is the change of free energy in 
the transition from gas to condensed phase 



^{kT/n) \n{P/Psat 



(4) 



where P and Psat are the partial pressure and saturation 
pressure of carbon. 

The evolution of partial pressure of carbon and gas tem- 
perature for models A2, A3, and B is shown in Fig. [5^. For 
the model Al, which is not shown, the behavior of calcu- 
lated values is intermediate between A2 and B. Using a stan- 
dard expression for Zeldovich factor (cf. Hasegawa & Kozasa 
1988) and the partial pressure of carbon from the interaction 
models one can compute from equations (2-4) the nucleation 
time (Fig. [SJj) for the adopted grain final radius a — 10~^ 
cm. It should be stressed that the result is not sensitive to 
adopted a. For a = 10~^ cm the value tnuc = 10'* s is at- 
tained by 0.8 day later. Apart from the standard mixture, 
we explored also the case of five times higher carbon and 
oxygen abundance He : C : O = 0.5 : 0.3 : 0.2 (Fig. [5j;,d). 
In this case the nucleation occures a bit earlier. 

A steep time dependence of tnuc{t) function permits 
us to recover the dust formation epoch for each particular 
model. Defining the dust formation epoch by the condition 
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Figure 5. Left panels: Loci of the radiation temperature and 
pressure in the CDS with the standard composition (upper panel) 
and five times larger carbon abundance (lower panel) . Thick solid 
line is for model B, intermediate line is for model A2 and thin line 
is for model A3. Tick marks correspond to every 5 days; the tick 
on day 50 is indicated. The phase equilibrium for graphite (dot- 
ted line) is taken according to Haines & Tsai (2002). Right panel: 
Evolution of nucleation time in the CDS with the standard com- 
position (upper panel) and five times larger carbon abundance 
(lower panel). Cases of the gas temperature equal to the radi- 
ation and the dust temperature are shown as solid and dashed 
lines respectively. The line thickness increases along the sequence 
A3, A2, B. 



tnuc = 1 day we find for the standard abundance that dust 
forms between days 42 and 54 for the low energy models and 
at 30-36 d for the high energy model A2. In the case of five 
times higher carbon abundance the dust forms at 40-48 d for 
the low energy models and at 28-32 d for the high energy 
model. The earlier dust-forming phase in the high energy 
model is related to the higher CDS velocity, i.e., lower gas 
temperature at the same age. This modelling shows that the 
dust forms in the CDS with the standard abundances of He, 
C, and O at the epoch close to day 50. 

The maximal amount of dust for the standard abun- 
dance is ~ 0.005 Mq assuming that all the oxygen is bound 
in CO molecules and the remaining carbon is converted into 
the dust. Around day 100 the CDS with the radius of ~ 10^*^ 
cm (Fig. 13)) and the dust mass of ~ 0.005 Mq has the opti- 
cal depth T ~ 500 in R band assuming extinction efficiency 
Qe = 2'jTa/X. The found upper limit is by a factor of ~ 10"^ 
larger than observational estimates (Mattila et al. 2008) 
which indicates either low efficiency of the dust formation, 
or clumpy structure of the dusty shell. Another interesting 
point is that the one-zone model considered above predicts 
rapid dust condensation, which is not the case: observations 
show slow increase of the amount of dust in SN 2006jc (Mat- 
tila et al. 2008). This means that conditions in the CDS 
matter perhaps are not uniform. A spacial variation of the 
gas temperature is probably a primary factor responsible for 
the gradual process of the dust condensation. 

We considered above the dust formation only in the 
CDS created by the reverse shock. However, one cannot rule 
out that the cool gas of shocked CS clouds in the forward 
shock might be also an appropriate site of the dust formation 
as well. With some modification this possibility corresponds 
to the scenario of the dust formation in the cool shell of the 
forward shock discussed by Mattila et al. (2008). 



5 DUSTY SHELL AND BLUE SHIFT 

The scenario proposed here suggests that He I emission lines 
originate from fragmented shocked CS clouds in the forward 
shock. To simulate line profiles affected by the dust absorp- 
tion one need to specify the distribution of emissivity and 
velocity of the line-emitting gas. I suggest that the unab- 
sorbed line profile is determined by the flow of cool frag- 
ments of shocked CS clouds in the range Ri < r < R2, 
where iii = 1 is the radius of the CDS, or more precisely, of 
the contact surface between shocked ejecta and shocked CS 
gas, and R2 is the radius of the forward shock front adopted 
to be 1.25 (cf. Chevalier 1982b). A possible collision of a 
cloud with the CDS would result in the sudden acceleration 
of the cloud shock; the shock gets adiabatic and cloud is 
demolished instantly. The resulting fragments in this case 
get hot and do not contribute in He I emission. The inte- 
grated emissivity is assumed to be constant in the range 
Ri < r < R2, which is sensible approximation, if the bulk of 
CS clouds crosses the postshock layer before they get com- 
pletely fragmented. 

The velocity distribution of the line-emitting gas gen- 
erated by CS clouds interaction with the forward shock is 
not trivial. The shocked cloud acquires a cometary structure 
with the low velocity core and high velocity tail of fragments 
(Klein et al. 1994). For particular cloud the flow kinemat- 
ics is approximately monotonic: the velocity increases from 
the core towards farmost tail. The smallest fragments have 
maximal velocities but they also are prone to rapid mixing 
with the hot gas (Klein et al. 1994). We assume that the 
smallest fragments of cool gas dissapear before they acquire 
the velocity of the forward postshock flow. From the mass 
conservation it follows that the average density of cool gas in 
the tail monotonically decreases with velocity. The combined 
flow of clouds at different stages of interaction cannot be re- 
duced to a single- valued function v{r). Instead, at the given 
radius a broad spectrum of velocities f{v) is expected with 
the minimal velocity corresponding to cloud shock speed 
vi ~ 10'^ km s~^ and maximal velocity Um(r) monotonically 
rising with r and determined by "oldest" shocked clouds 
that are colliding with the CDS. I adopt the linear law for 
the maximal velocity Vm — vi + {v2 — vi)(r — Ri), where 
V2 is a fitting parameter. The velocity density distribution 
function at the given radius is taken in the simple form 

/(„)^C(^^^), (5) 

\V2 — Vl ./ 

where C is normalizing factor. 

The CDS is subject to Rayleigh- Taylor instability, frag- 
mentation, and mixing in the forward shock layer (Chevalier 
1982a; Blondin & Ellison 2001), so the dusty CDS material 
should be spread between the contact surface Rd^i = Ri — 1 
and some external radius Rd^2 which is a free parameter. The 
emission line profiles affected by the dust is simulated here 
using Monte-Carlo technique. The dust in our model has ve- 
locity of the CDS for which we adopt 9000 km s"'^. This pa- 
rameter weakly affects the scattered radiation. The Rayleigh 
phase function and scattering albedo u = 0.6 are assumed. 
I checked different values of u and found only minor dif- 
ferences which refer to the intensity of the broad scattered 
component. Note, this component, although weak, may be 
revealed in the polarized light, if the asymmetry of dust dis- 
tribution is significant. I also tried the phase function dom- 



Circumstellar interaction in type Ihn supernovae and SN 2006jc 7 



Table 2. Parameters of the line profile model 
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inated by back scatter which would be the case, if the dust 
resides in optically thick clumps. The difference, primarily 
in scattered component, turns out small. 

The optimal model profiles are plotted (Fig. [6]) together 
with observed He I 7065 A line profile of SN 2006jc for dif- 
ferent epochs (Foley at al. 2007). The fit is good for the 
three first dates and sensible in the last plot. This can be 
considered as success of the proposed model. The recovered 
parameter V2, extinction optical depth and radii of the dusty 
shell are given in Table 2. The decrease of V2 suggests that 
the contribution of the high velocity line-emitting gas falls 
with time. Note, the parameter V2 is determined by unab- 
sorbed blue part of profile and does not depend on the dust 
distribution and the optical depth. It is not clear whether 
the decrease of the high-velocity component is related with 
the cloud fragmentation and fragment acceleration or with 
the mixing of high velocity fragments in the hot gas. There 
is a hint that the relative radius of the outer boundary 
of the dusty shell Rd,2 increases with time. Interestingly, 
such a behavior qualitatively agrees with the initial evolu- 
tion of a mixing layer driven by the Rayleigh- Taylor insta- 
bility (Chevalier & Blondin 1995; Blondin & Ellison 2001). 
A detailed comparison with published numerical models is, 
however, not possible, because the models for SN 2006jc pre- 
sented here do not have power-law density profiles for the 
ejecta and the wind. 

I explored cases with larger inner radius of the dusty 
shell, Rd,i > 1, but found that the fit worsens in this case. 
This means that in the preferred model the center of the 
dust layer is shifted inward relative to the center of the line- 
emitting layer. This is expected in the proposed picture in 
which most dust originates from the fragmented CDS and 
the line-emitting gas is related with the shocked CS clouds. 



6 DISCUSSION 

The proposed scenario of SN Ibn with the strong early CS 
interaction turns out successful in explaining some unusual 
observational properties of SNe Ibn. However, until now we 
ignored some questions that might have interesting conse- 
quences for our model. It is worthy, therefore, to pose and 
briefly discuss challenging issues and, if possible, to impose 
additional constraints for the model. 



6.1 What does power blue continuum? 

The proposed interpretation of the blue bumpy continuum 
suggests that the most emergent optical radiation originates 
from the CDS. In case, when the interaction luminosity dom- 
inates, the energy from the reverse and forward shocks is 
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Figure 6. Line profile of He I 7064 A in SN 2006jc spectra for 
four epochs. The model and observations (Foley et al. 2008) are 
shown as a solid line and dotted line, respectively. 



deposited into the CDS in a usual way via absorption of X- 
rays and, perhaps, the thermal conductivity in the mixing 
zone of the forward shock. 

Another possibility could be that the dominant source 
is ^^Ni decay. In this case the energy is deposited into CDS 
owing to the Compton scattering of gamma-quanta and the 
absorption of ultraviolet radiation emitted by the unshocked 
ejecta in turn powered by the radioactive decay. The effi- 
ciency of the gamma-ray deposition to the CDS is of the 
order of the ratio of the CDS mass to the mass of ejecta, 
which is small unless ejecta mass is low < 1 Mq. 

Alternatively, the CDS gas could be ionized and excited 
by the absorption of the ultraviolet radiation from the un- 
shocked ejecta. The problem with this mechanism is that in 
the unshocked freely expanding SN Ibc envelope the ultravi- 
olet radiation is strongly reprocessed by the photon scatter- 
ing and splitting in numerous ultraviolet metal lines, so the 
emergent ultraviolet fiux gets strongly suppressed relative to 
visual. The situation in this respect is similar to SNe la (cf. 
Pinto $ Eastman 2000). At first glance the ultraviolet fiux 
might become stronger, if the ejecta mass is low. However, 
for this to be the case the mass should be much lower than 
that of SNe la, i.e., M < 1 Mq. We thus conclude that in 
the radioactive model the CDS can be the source of the blue 
continuum, if the ejecta mass is low, < 1 Mq . 



6.2 Where does He I absorption form? 

Helium lines (e.g. 3889 A) in the blue part of SN 2006jc spec- 
trum show P Cygni profiles with the absorption minimum 
at ~ —1000 km s~^ and a comparable equivalent width of 
absorption and emission components (Foley at al. 2007). A 
natural explanation of this phenomenon is that the contin- 
uum in blue region is strong, so the scattering component 
dominate or comparable with the net line emission. In fact 
absorption components are also seen in O I 7773 A and m 
He I 6678 A lines (Anupama et al. 2008). The maximal ve- 
locity observed in the blue absorption wing of He I 3889 A 
line is ~ 1900 km s~^, three times larger than the wind 
velocity, v ~ 600 km s~^. Apart from the wind, therefore. 
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an absorbing high velocity component (> 600 km s^^) is 
needed to account for the line absorption. 

Three possibilities are conceivable for the line- absorbing 
gas in the range of 600-2000 km s"^: (a) the CSM in fact has 
the free expansion kinematics v oc r with the most material 
in the range of 600-2000 km s"'^; (b) the hne-absorbing gas 
is simply the same line-emitting gas in the forward shock; 
(c) the CSM is flow with the constant velocity of ~ 600 km 
s~^ but in the preshock zone the CS gas is accelerated up 
to 1900 km s^^. The first option predicts that the emission 
lines related with the CSM should get broader with time. 
In fact, the absorption of He I 3889 A does not show this 
trend according to spectra taken by Foley et al. (2007). The 
second option suggests that fragments of CS clouds in the 
forward shock should have large covering factor ~ 1, which 
is unlikely, although is not rulled out. The third possibility 
is similar to that suggested for the Ha high velocity CS ab- 
sorption component in early SN 1998S (Chugai et al. 2001), 
where this feature was attributed to the preshock gas accel- 
erated by either radiation or cosmic rays generated in the 
forward shock. At the distance of ~ lO'^^ cm the radiative 
force in case of SN 2006jc is weak and cannot accelerate the 
preshock gas up to more than 100 km s"'^ in the wind frame 
and this mechanism should be rejected. 

As to cosmic rays, the situation is more optimistic. The 
point is that in case of clumpy CSM the cosmic ray pressure 
is build up by the overall shock kinetic luminosity, while 
the acceleration is applied primarily only to the intcrcloud 
gas with a small mass fraction of the inflowing CSM. The 
intercloud gas, therefore, can acquire significamt velocity in 
the preshock zone. Let pic is the density of the intercloud 
gEis that is substantially lower than the average CS density 
p. Momentum and mass conservation result in the maximal 
preshock velocity of the intercloud gas 

u=\4>Vs—, (6) 
6 Pic 

where i;, is the forward shock speed, 4> is the ratio of the 
energy density of cosmic rays to the density of the shock 
kinetic energy 0.5pwf . For Vs = 9000 km s~^ and p/pic = 5 
one needs <j} ~ 0.2 to produce the preshock velocity of tt « 
1300 km s^^ in the wind frame. This estimate demonstrates 
that velocities of absorbing gas up to 1900 km ~^ might be 
explained, if the line- absorbing gas is associated with the 
intercloud gas accelerated in the cosmic ray precursor. 



6.3 Why intensities of Mgll and OI lines are 
similar? 

Between days 22 and 60 intensities of Mg II 7889 A and O I 
7773 A lines in SN 2006jc are comparable (cf. Anupama et 
al. 2008). Given similar ionization and excitation potentials 
of these species a straightforward implication of this fact 
is that the abundances of oxygen and magnesium are com- 
parable. This indicates that we perhaps sec the matter in 
which oxygen is strongly depleted by CNO burning. Indeed, 
massive stars at the boundary between helium core and hy- 
drogen envelope have a zone, in which oxygen is depleted 
owing to CNO burning by a factor of 10-20 (cf. Hirschi et 
al. 2004; Meynet & Maeder 2003). The low abundance of 
hydrogen in the CSM around SN 2006jc is consistent with 



the conjecture that we might see the matter with strongly 
depleted oxygen. 

An alternative explanation is that Mg/0 ratio is cosmic, 
and intensities of Mg II and O I lines are comparable because 
these lines are thermalized, i.e., level populations are close 
to Boltzmann distribution and lines are optically thick. The 
thermalization in a certain line occurs if e2iNs > 1, where 
Ns is the average number of conservative scattering in the 
line before the photon escape, £21 = <l2i'rie/A2i is the ther- 
malization parameter, and q2i is the coUisional de-excitation 
rate. Using estimates based on the CSM density, pressure 
equilibrium. He I line intensity, and Boltzmann excitation of 
the lower level I estimate ^21 Ns ~ 10 for the Mgll 7889 A 
line on day 22. This confirms the possibility of the thermal- 
ization of Mgll and OI emission lines. 

What explanation of the comparable intensities of Mg II 
and 01 lines — similar abundances of oxygen and magne- 
sium or thermalization — is correct remains an open ques- 
tion. The conjecture of the similar abundances of Mg and O 
seems to contradict to the fact that the narrow absorption 
is absent in Mg II 7889 A, while it is present in O I 7773 A 
(cf. Anupama et al. 2008). 

6.4 Ejecta, presupernova and explosion 

Our modelling of light curve and other properties of 
SN 2006jc does not permit us to constrain the ejecta mass 
and explosion energy, although models with the moderate 
energy (1 — 2) x 10®^ erg are preferred, if one admits that 
the dust forms in the CDS with the velocity of ~ 10* km 
s"^ Moreover, models of SN 2006jc with and without '^^Ni 
are plausible, although both have its own drawbacks. The 
model of a standard SN Ibc powered by ^^Ni decay faces 
problems if we identify the blue continuum with the emis- 
sion of the CDS. Only the low mass ejecta < 1 Mq could 
reconcile the radioactive model with the proposed scenario 
of the blue continuum. On the other hand the drawback 
of the model without ^^Ni is that it predicts unacceptably 
strong X-ray luminosity. 

A more certain judgement on the source of SNe Ibn 
luminosity at the nebular epoch could be made in principle 
on the basis of the late time (t > 200 d) photometry and 
spectra of SN Ibn. A deviation of the bolometric light curve 
from that predicted by the radioactive model would favour 
the CS interaction as a primary energy source. On the other 
hand uniformity of light curves and their correspondance 
to that expected for the radioactive model would favour the 
radioactivity. If the radioactivity as a primary energy source 
will be confirmed the second crucial question would remain, 
whether the ejecta is a standard SN Ibc or the object similar 
to our low mass model A3 with ~ 1 Mq ejecta and '"''Ni mass 
of ~ 0.3 Mq. This issue could be elucidated by the late time 
spectra that may or may not reveal strong [O I] 6300, 6364 
A emission, the canonical feature of SN Ibc, at the nebular 
stage. The strength of this emission relative to Fell line 
continuum would confirm or reject a picture of explosion of 
a standard SN Ibc slighly spoiled by CS intraction. The dust 
absorption could hamper this test, although it may well be 
that due to dumpiness of the dusty shell the attenuation 
will not be large at the late nebular epoch. 

The CSM mass around SN 2006jc at the distance r < 
2 X 10^^ cm in our models (Table 1) is in the range of 
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0.02 - 0.05 Mq. With the outflow velocity = 600 km 
s^^ this mass had to be lost by pre-SN during the last year 
before the explosion with the rate of ~ 0.02 — 0.05 M© yr~^ 
and kinetic luminosity of ~ (2.5 — 6) x 10^® erg . The out- 
flow characteristics suggest extraordinary behavior of prc- 
SN. The pre-SN mass loss might be related with the bright 
outburst two years before the SN 2006jc explosion (Nakano 
et al. 2006; Foley et al. 2007; Pastorello et al. 2007). How- 
ever, if the ejection episode was as brief as the optical out- 
burst then the CS shell should be narrow and reside at the 
distance of about 4 x 10^^ cm in contrast to our model which 
suggests a significant fraction of CSM in the close vicinity, 
at r < 2 X 10^'' cm. Assuming the CS shell expansion veloc- 
ity of 2400 km s^^ one gets the CS shell oven farther out, 
at the distance of ~ 1.6 x lO" cm (Mattila et al. 2008). 

A problem of a vigorous mass loss that preceeds the SN 
explosion is familiar in relation with SNe Iln. Following the 
conjecture invoked for SN 1994W (Chugai et al. 2004b) one 
might suggest that in some cases, including SNe Ibn, the 
mass ejection by prc-SN is initiated by the flash of nuclear 
burning, e.g., of neon. However, the current theory of stellar 
evolution does not predict the violent mass loss of pre-SNe 
with initial masses < 100 A^q (Hcgcr et al. 2003). Only mod- 
els of massive stars in the range of 100 — 140 M© reveal the 
instability that could result in the strong nuclear outbursts 
and mass ejection (Heger et al. 2003). The physics is in pair 
production, which leads to infall, explosion of carbon and 
oxygen followed by mass ejection. These stars can experience 
several ejection episodes before the ultimate collapse and 
collisions between consecutively ejected shells could produce 
optical outburst resemblant supernovae events (Heger et al. 
2003). Recently, extremely luminous SN 2006gy was identi- 
fied with a similar phenomenon (Woosley et al. 2007). It is 
not clear, whether such a scenario is applicable to SNe Ibn. 
If does, then these supernovae should be attributed to the 
CS interaction in the absence of "''''Ni. An important impli- 
cation of this scenario is that the luminous massive helium 
star should remain at the position of SN Ibn. 

An alternative scenario, in which SN explosion might 
be preceeded by the vigorous mass loss, is a binary system 
in which a black hole experiences a merging with the helium 
core of a red supergiant. The supernova explosion in this case 
is presumably driven by the non-relativistic jets generated 
by the rapid disk accretion of the helium core material onto 
the black hole. In fact this possibility is a version of a merger 
scenario of gamma-ray bursts (Fryer & Woosley 1998; Fryer 
et al. 1999). The merging of a neutron star with a red su- 
pergiant leads to the similar scenario because the inspiraling 
neutron star that accrete at the Bondi-Hoyle accretion rate 
grows black hole already in the hydrogen envelope (Cheva- 
lier 1993; Fryer et al. 1999). The binary scenario predicts 
asphericity of ejecta caused by jets and of the CSM caused 
by the equatorial outflow from the common envelope. 



7 CONCLUSIONS 

The goal of this study was to present a model that might 
account for the unusual properties of SNcIbn including 
SN 2006jc. I argue that the CS interaction with a dense 
CSM is essential energy source for the early light curve of 
SN Ibn. I demonstrated that the CS interaction at the early 



stage loads to the formation of the CDS which in turn is re- 
sponsible for the smooth continuum of early SN 2000er and 
unusual blue bumpy continuum of SN 2006jc. I showed that 
the dust can form in the CDS at about day 50 in accord with 
observations. The modelling of the He I line profile affected 
by the dust absorption confirms a picture in which the line- 
emitting gas is associated with the CS clouds shocked and 
fragmented in the forward shock, while the dust is related 
with the fragmented CDS. 

It should be stressed that scenario in which SN Ibn su- 
pernova is the explosion of a standard or luminous SN Ibc 
sligtly modified by the CS interaction remains unconfirmed 
and photometry and sectroscopy of SNe Ibn at the late neb- 
ular stage are needed to clarify the issue. 
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